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The collapse and swelling behavior of a generic homopolymer is studied using implicit-solvent, 
explicit- cosolvent Langevin dynamics computer simulations for varying interaction strengths. The 
systematic investigation reveals that polymer swelling is maximal if both monomer- monomer and 
monomer-cosolute interactions are weakly attractive. In the most swollen state the cosolute density 
inside the coil is remarkably bulk-like and homogenous. Highly attractive monomer-cosolute inter- 
actions, however, are found to induce a collapse of the chain which, in contrast to the collapsed case 
induced by purely repulsive cosolvents, exhibits a considerably enhanced cosolute density within the 
globule. Thus, collapsed states, although appearing similar on a first glance, may result from very 
different mechanisms with distinct final structural and thermodynamics properties. Two theoretical 
models, one based on an effective one-component description where the cosolutes have been inte- 
grated out, and a fully two-component Flory - de Gennes like model, support the simulation findings 
above and serve for interpretation. In particular, the picture is supported that collapse in highly 
attractive cosolvents is driven by crosslinking-like bridging effects, while the ratio of attraction width 
to cosolute size plays a critical role behind this mechanism. Only if polymer-cosolute interactions 
are not too short-ranged swelling effects should be observable. Our findings may be important for 
the interpretation of the effects of cosolutes on polymer and protein conformational structure, in 
particular for highly attractive interaction combinations, such as provided by urea, GdmCl, Nal, or 
NaC104 near peptide-like moieties. 



I. INTRODUCTION 

Single polymer coils in solution swell in good solvents 
and collapse in bad solvents^ In a good solvent the effec- 
tive interaction between polymer monomers is repulsive 
which tends to swell the polymer. In a bad solvent the 
interaction is essentially attractive and the coil shrinks 
until hindered by steric packing effects. On the simplest 
theoretical level these trends can be qualified by effec- 
tively one-component mean-field treatments as pioneered 
by Flory and de GenneSfi^ Here the system free energy 
F(R) in dependence of coil size R is typically written 
as2£ 

F(R) ~ R 2 /Nb 2 + B 2 N 2 /R 3 + B 3 N 3 /R 6 , (1) 

where the first term is the elastic energy of an ideal chain 
with N monomers and segment length b, and the other 
terms account for mutual monomer interactions globally 
expressed by an virial expansion up to second order in 
density ~ N/R 3 . The virial coefficients have to be con- 
sidered as based on effective pair potentials because the 
(co)solvent degrees of freedom have been integrated out. 
Scaling law predictions from (1) are basis for the discus- 
sion of collapse and swelling in polymer science^^ 

The problem of such a perspective is that the infor- 
mation about the detailed effects of solvent or cosol- 
vents are lost. In the recent years, however, the com- 
plexity and resolution of experimental investigations has 
increased, and there is growing interest about how the 
specific binding mechanisms of cosolutes may alter global 
polymeric properties, such as the coil-globule transition. 4 
A particularly important example is the specific effect 



of cosolutes, such as ions, osmolytes, or denaturants on 
protein structure and stability. Here, general poly- 
mer principles are typically employed as a starting point 
for the investigation of those cosolute effectsi 1 ^! In this 
spirit, relatively simple thermosensitive polymers such as 
the poly(N-isopropylacrylamide) (PNIPAM) homopoly- 
mer often serve as a experimental peptide model . 18 ' 19 
On the computational side, simulations allow to ad- 
dress fundamental questions in the framework of poly- 
mer physics, such as swelling and collapse of idealized, 
purely hydrophobic homopolymer in denaturants, os- 
molytes, or saltS)22r— or more generic polymer-cosolute 
systems £2r— 

An important experimental reference to our study are 
the works by the Cremer group on (methylated) urea and 
ion-specific effects on the lower critical solution temper- 
ature (LCST) of PNIPAMi&±2 and the relatively sim- 
ple, elastin-like peptides! 19 ' 27 The qualitative change in 
the LCST, AT(p), with cosolute concentration p can be 
used as an index whether polymer coils swell or collapse 
with the addition of cosolutes^ For instance, PNIPAM 
swells in methylated urea, while it collapses with the ad- 
dition of nonmethylated urea. In contrast, the relatively 
hydrophobic elastin-like peptides swell in urea, consistent 
with computer simulations of simple homopolymersi22r— 
For both PNIPAM and elastin peptides, both being clec- 
troneutral, collapse and swelling in salt strongly depends 
on salt type. Hence, there is a wide variety of observed 
behaviors which are challenging to categorize. 

Typically polymer collapse (or protein stabilization) 
is argued to originate from the preferential exclusion of 
repulsively interacting cosolutes, thereby trying to mini- 
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mize cosolute- accessible surface area^ i 10 ' 11 ' 28 Sagle et at, 
however, argue that urea collapses PNIPAM due to a di- 
rect binding mechanism featuring strong, short-ranged 
attractions provided by multivalent hydrogen bonds. 19 
Crosslinking then leads to shrinking of the coil. This 
were not the case, they say, for methylated urea, or urea 
interacting with peptides, where interactions are only 
monovalent and very weakly attractive. Hence, specific 
values of binding parameters such as attraction strength 
and width are probably decisive whether a given polymer 
swells or not. 

The crosslinking effect may also be important for other 
cosolutes or salts, such as NaClC>4. Its effects on the 
LCST of PNIPAM does not conform with the usual 
Hofmeister seriesi*^ 8 - On the other hand, it is known to 
strongly interact with the peptide group . 18 ' 19 Indeed re- 
cent explicit-water computer simulations together with 
circular dichroism (CD) and Forster resonance energy 
transfer (FRET) measurements on NaC104-destabilized 
a-helical peptides demonstrated highly compact, disor- 
dered states crosslinked by a network of sodium and 
perchlorate ions^ 9 - Thus, collapsed states, although ap- 
pearing similar on a first glance, may result from very 
different mechanisms with quite distinct final structural 
and thermodynamics properties. Indications of the latter 
stem from the experimental characterization of NaC104 
denatured molten globule o 5 ' 29 ' 39 , which may be cosolute 
rich, with implications in protein folding . 31 ' 32 We note 
that similar considerations may be important in the ion- 
induced collapse of polyelectrolytes by ion condensation 
beyond simple Debye-Huckel electrostatics.— ~— Hence, 
finding minimalistic models which can describe this vast 
variety of effects are in need. 

The aim of this paper is to theoretically investigate 
polymer collapse and swelling under the action of coso- 
lutes on a highly generic level. For that we first start to 
employ implicit-solvent computer simulations of a Gaus- 
sian polymer chain including explicit cosolutes, and we 
systematically vary monomer-monomer and monomer- 
cosolute attraction strengths. Our polymer model is sim- 
ilar to that of Toan et at— who systematically investi- 
gated polymer swelling and collapse (with consequences 
to conformational kinetics) in dependence of monomer- 
monomer attraction strength in purely implicit solvent. 
Our simulations indicate chain collapse for highly at- 
tractive cosolute conditions as also found in recent on- 
and off-lattice computer simulations by Antypov and El- 
liott The results from the simulations are then ratio- 
nalized by two theoretical models: First, an effective 
interaction model where the action of the monomer- 
cosolute is explicitly integrated out in a statistical me- 
chanics framework^ The results are then discussed on 
a one-component level as in eq (I). Secondly, we ex- 
tend the Flory-deGennes description in eq. (1) to a full 
2-component description and calculate swelling behav- 
ior in dependence of interaction strengths and cosolute 
density. Both theoretical descriptions qualitatively agree 
with the simulations giving important rational. Impor- 



tantly, all three approaches point to polymer collapse at 
highly attractive cosolute conditions due to crosslinking- 
like bridging effects as argued by Sagle et a/.— 

II. POLYMER-COSOLUTE LANGEVIN 
COMPUTER SIMULATION 

A. Model and Methods 

In our simulation we consider a single, coarse-grained 
homopolymer with N = 100 connected monomers in an 
implicit solvent background in a volume V. In addition, 
N c explicit cosolutes with mean number density p 9 = 
N c /V are added to the system. All particles, polymer 
monomers and cosolutes, are interacting via the Lennard- 
Jones (LJ) pair potential 

V ij (r)=Ae ij [(a ij /r) 12 -(a ij /r)% (2) 

where i = m,c stands for monomer or cosolute atom, re- 
spectively. The value of the size <xy = a = 0.3385 nm is 
chosen to be be fixed for all interactions in our simula- 
tions. It is the same value as in Toan et aZ.'s work— and 
is comparable to the typical size of the methyl group. 
The interaction energy will be systematically var- 
ied as described below. The polymer is modeled by a 
Gaussian chain with harmonic nearest-neighbor bond in- 
teractions Vbond.ij = 0.5fe(rj — Vj) 2 and a spring con- 
stant k = 320 kJ mol -1 nm~ 2 . The latter is chosen 
such that we end up with the same effective Kuhn length 
b = 0.38 nm for the ideal polymer as in Toan et aVs 
work.' 5 ' The mean end-to-end distance in the ideal case 
is thus R — \f~Nb — 3.8 nm. The corresponding ideal 
radius of gyration i? 9 = 1.55 nm. The LJ interaction 
between a monomer and its next two nearest neighbors 
is excluded. 

The polymer-cosolute system is simulated in the NVT- 
cnscmble using stochastic computer simulations. Every 
atom is propagated via the Langevin equation 

mr, + mtti = ^ F y + , (3) 
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where rj is the position of a particle i, m = 8 amu 
its (irrelevant) mass, £ = 0.5 ps _1 is the friction con- 
stant, Fij the force between particles i and j, and 
the stochastic force stemming from the solvent kicks. 
The stochastic force has zero mean (F^) = and 
fluctuation-dissipation is obeyed via (F^ R ' (t)Fj\t')) — 
2kBT£,5ij6(t — t'). The simulations are performed using 
the Gromacs simulation package^ 9 , with periodic bound- 
ary conditions, an integration time step of 4 fs, and lines 
bond constraints. The simulation box is cubic and has a 
box length of 13 nm. We fix the mean cosolute density 
p 9 = 3 nm~ 3 (~ 5 mol/1) chosen to represent a typical 
denaturant density. With that we end up with N c = 6591 
cosolutes in the simulation. We simulate every system for 
at least 1 /its. 



If the polymer is in a swollen or collapsed state can be 
judged by inspection of the radius of gyration R g scaled 
by the one of the ideal polymer R°. If R g /R g > 1 the 
polymer is swollen, if R g /R g < 1, it is collapsed. The 
statistical error of R g in our simulations is estimated by 
block averages: the trajectory is divided into n blocks of 
the same length (at least 200ns, but up to 1/xs) and the 
standard deviation of these block averages with respect 
to the mean was calculated and divided by ^Jn to obtain 
the error. 

Whether the LJ interaction Vjj(r) is globally attrac- 
tive or repulsive can be quantified by the second virial 
coefficient defined via 

B 2=~\ J d 3 r[cxp(-/3^(r))-l], (4) 

where j3 = (fc^T) -1 is the inverse thermal energy and 
we use fcgT = 1 as the energy scale in the following. 
For the LJ potential, the B 2 in dependence of the LJ 
interaction energy e is shown in Fig. 1: for interaction 
values e < 0.3, the B 2 is larger zero since excluded vol- 
ume contributions from distances r < a dominate. At 
e = 0.3 the B 2 vanishes, i.e., interactions are ideal on a 
2-body level. For e > 0.3 effectively the interactions are 
attractive. In the polymer-only case (no cosolutes), those 
regimes can be identified with the usual good solvent 
regime (swelling), 8-solvent (ideal), and bad solvent (col- 
lapse), respectively— Indeed previous simulations showed 
a collapse transition of the lOOmer in the weakly attrac- 
tive regime with e around 0.4. For e > 0.6 the chain 
collapses in dense, compact states. This regime was char- 
acterized as highly attractive^ 

In our simulations we vary systematically both e mm 
and e mc between and 1, thereby studying the whole 
range of repulsion to strong attraction for both compo- 
nents, that is, monomer- monomer and monomer-cosolute 
interactions. For simplicity the cosolute-cosolute inter- 
action will be fixed to e cc = 0.3 mimicking a near-ideal 
cosolutc solution. Nonideality may also influence poly- 
mer swelling and shrinking behavior due to osmotic ef- 
fects but these effects are out of scope of this paper. 
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FIG. 1: Second virial coefficient eq. (J4j) of the Lennard- Jones 
interaction eq. (J2J) versus LJ interaction energy e. The B2 
value can be used to classify the interaction into mainly 're- 
pulsive' (B2 > 0), 'attractive' (B2 < 0), and 'highly attrac- 
tive' (B 2 /a 3 < 2). 
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FIG. 2: The squared radius of gyration R 2 , of the polymer 
scaled by the ideal value (Rg) 2 versus the monomer-cosolute 
interaction strength e mc for varying monomer-monomer in- 
teraction e mm . If (R g /R°g) 2 > 1 the polymer is swollen, if 
(Rg/R^) 2 < 1, it is collapsed. Error bars are typically of 
symbol size. 



B. Simulation Results 

In Fig. 2 we plot the size of the polymer expressed 
by the radius of squared gyration R g versus monomer- 
cosolute attraction e mc for varying monomer-cosolute at- 
traction e mm . Note that R g is scaled by (R g ) 2 , the size of 
the ideal polymer, to easily distinguish between swollen 
(R g /R g > 1) and collapsed (R g /R g < 1) states. Let us 
first focus on vanishing e mc = 0, i.e., the polymer-only 
case: the polymer is swollen in a good solvent e mm < 0.3, 
almost ideal for e mm — 0.4, and collapsed for e mm > 0.5. 



If now the cosolutes are 'switched-on' (e mc > 0) to 
a value of e mc = 0.1, where monomer-cosolute interac- 
tions are repulsive (cf. Fig. I), the chain shrinks con- 
siderably for e mm < 0.6. This is expected since the ex- 
clusion of repulsive monomers from the polymer region 
is entropically disfavored and leads to chain shrinking. 
For larger e mm > 0.6 the effects are small. For an in- 
creased e mc = 0.2 the polymer swells for all of the values 
of e mm . Hence, the solvent quality overall improves with 
less monomer-cosolute repulsion (cf. Fig. 1). For very 
strong monomer-monomer attractions, e mm > 0.6, the 
effect is quite small and the polymer stays essentially in 
the collapsed state. 

For further increasing e mc > 0.3, roughly where 
monomer-cosolute repulsion turns into attraction, the 



swelling of the polymer continues for all e mm to a max- 
imum value BJ^ ax whose corresponding e™® x value de- 
pends on e mm . For larger monomer- monomer attrac- 
tion e™® x is shifted to larger values. For values larger 
than e™° E the chain collapses again, i.e., there is a reen- 
trant collapse transition for an increasingly attractive 
monomer-cosolute interaction. This continuous crossover 
from a good solvent to a bad one with increasing at- 
traction has been observed already in a previous simu- 
lations.— For not highly attractive intrapolymer interac- 
tions e mm < 0.6 and highly attractive polymer-cosolute 
interactions e mc > 0.6, the chain can be even more 
strongly collapsed than in the case of highly repulsive 
monomer-cosolute (e mc = 0.1) interactions. Hence, this 
regime shows strong compaction by a highly attractive 
cosolute. 

The cosolute-induced effect is strongest for e mm ~ 0.4, 
that is, in the ideal to weakly attractive polymer regime. 
Overall the polymer crosses from collapsed (e mc ~ 0.1) to 
swollen (e mc ~ 0.3 — 0.8) to again collapsed (e mc > 0.8) 
states with the most dominant changes when compared 
to the other monomer- monomer interactions. Simulation 
snapshots of the polymer-cosolute system for e mm ~ 0.4 
and e mc =0.1 (collapsed), 0.5 (swollen), and 1.0 (col- 
lapsed) are shown in Fig. 3. Note the considerably dif- 
ferent amount of cosolutes in the vicinity of the two col- 
lapsed states (a) and (c). 




FIG. 3: Simulation snapshots of the polymer-cosolute system 
for e mm ~ 0.4 and (a) e mc =0.1 (collapsed), (b) 0.5 (swollen), 
and (c) 1.0 (collapsed). Cosolutes (yellow spheres) are shown 
which are found in a cut-off distance of ct2 1//6 = 3.8A to the 
polymer backbone (connected red spheres). 

To further characterize structural details in the regimes 
described above we plot the radial one-particle density 
distributions of the monomers and cosolutes Pm{r) and 
p c (r) in Fig. 4. With r we denote the radial distance 
to the center-of-mass of the polymer chain. The profiles 
are plotted for a fixed e mm = 0.4 and varying e mc = 0.1, 
0.5, and 1.0. In the collapsed states the polymer density 
profiles are similar. However, the cosolutes are strongly 
depleted in the repulsive case, e mc =0.1, while their pres- 
ence is massively enhanced in the highly attractive case 
with e mc =1.0. At the most swollen polymer state for 
£mc = 0.5 the cosolute density is remarkably homoge- 
neous and bulk-like. 

III. EFFECTIVE INTERACTION MODEL 

In this section we attempt to rationalize some of the 
swelling and collapse trends we found above in the per- 
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FIG. 4: Radial density profiles pi(r) of polymer monomers 
(i = m) and cosolutes (i = c) around the center-of-mass of 
the polymer. Interaction values are e mm = 0.4 and e mc = 0.1, 
0.5, and 1.0 (see legend) corresponding to the snapshots in 
Fig. 3. 



spective of an effective one-component model. For this 
the effects of the cosolutes have to be integrated out to 
derive an effective cosolute-induced interaction between 
two monomers. This is reflected then in an effective B2 
coefficient as used in the one-component Flory approach 
eq (fTJ . For simplicity we will model the monomers as pla- 
nar plates and aim only at qualitative statements. We 
borrow thereby from a similar model introduced within 
the framework of mesosurface attraction induced by ad- 
hesive particles 

Consider two planar surfaces (monomers) with area A 
in a surface-to-surface distance d in contact with a reser- 
voir of cosolutes at concentration p° c . The cosolutes in- 
teract with the surface with the generic potential V(z) 
as shown in Fig. 5: the cosolutes are hard spheres with 
a diameter a and have an additional attractive (adsorp- 
tion) energy of strength e < and not too large width 
5 < a/2. Now consider three situations l)-3) as also 
depicted in Fig. 5. In situation 1) the plate distance is 
d < it, i.e, no cosolutes fit between them. This situa- 
tion can be coined 'depletion' situation. In situation 2), 
d > 2a, and cosolutes can solvate both of the surfaces. In 
situation 3), d ~ a, one bound particle can interact with 
the two surfaces simultaneously This situation can be 
named 'bridging' or 'crosslinked' situation. Interactions 
between the cosolutes, e.g., packing effects, are neglected. 
These effects have been included in a more complete cal- 
culation within a similar model previously^ 

In case 1) the proximity of the plates increases the ac- 
cessible volume for the cosolutes. The latter gain config- 
urational entropy proportional to the freed volume Aa. 
This leads to a favorable (grand canonical) free energy 
per area 

ni/A = -ap° c (5) 
for d < a. That is the well-known effective surface at- 
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FIG. 5: Top: Generic interaction potential V(z) between one 
cosolute and one plate-like monomer in a center-of-mass to 
surface distance z as assumed in our statistical mechanics ap- 
proach. a/2 is the cosolute hard-core radius, while e < is 
the attractive well depth and S is its width. Bottom: sketch 
of the three interaction situations of the simplified plate- 
like monomers (blue plates) in presence of the cosolutes (red 
spheres) . 



traction by particle depletion. 40 The grand canonical free 
energies for situations 2) and 3) can be explicitly written 
down and are 



and 



^2 = N h 
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(7) 



respectively, where A is the thermal wavelength, /z the 
chemical potential, and Nh is the number of bound par- 
ticles. Note the difference between (2) and (3): In case 
2), a particle gains energy e by binding to one surface, 
but 2 surfaces, that is, twice the configuration space as 
in 3) is available in total (hence the factor 2 before A in 
the log-term). In situation 3), one particle gains 2e by 
binding to two surfaces simultaneously, but effectively 
only one plane for the cosolute translation is available. 
The energies are minimized by the Boltzmann equations 
N/(2SA) = p°exp(-e) and N/(SA) = p°exp(-2e), re- 
spectively, where we used p Q c = A _3 exp(p). Plugging 



the solution back in (1) and (2) we obtain the minimum 
free energies per area 



and 



Q.* 2 / A = -25p° c exp(-e) 



fl* 3 /A = -<5p°exp(-2e). 



(8) 



(9) 



Inspection of the results shows that for very small 
|e| situation 1) is favored over 2), that is, the deple- 
tion interaction wins. For increasing e, however, situ- 
ation 2) is favored over 1); when comparing i7* = f2| 
we find e* 2 = — lna/(2S) which depends on the inter- 
action range 5. For a LJ interaction, as in our simula- 
tion model, the width of the minimum roughly S ~ a/3 
which yields e* 2 — — 0.4 for the crossover from depletion 
attraction to full solvation of the two monomers. Since 
the minimum free energies of the two situations 1) and 
3), where d < a, are both higher, this implies repulsion 
between close monomers. In colloidal physics this effect 
was coined 'repulsion through attraction' ] 41 ' 42 The free 
energy for situation 3), where bridging is favored, de- 
creases faster with e than in situation 2), because of the 
2 in the exponent. By equating (jHJ and (j9|), we obtain the 
crossover energy e| 3 = — In 2 ~ —0.7. For high attraction 
consequently situation 3) is favored. 

Thus, we can now generally identify three regimes 
depending on the attractive interaction strength e: i) 
|e| < 0.4: depletion attraction between surfaces (favor- 
ing collapse); ii) 0.4 < |e| < 0.7: repulsion through at- 
traction (favoring swelling); |e| > 0.7: polymer-polymer 
attraction by bridging (favoring collapse). The situation 
is summarized in the 'phase diagram' in Fig. 6 for the 
specific value of 8 = a/3. Note, however, that the at- 
traction width plays a nontrivial role: if the attraction 
range is small compared to the cosolute size, S < cr/4, 
then |e* 2 | > 0.7, and situation 2) becomes metastable. 
In contrast, if S is not small, 5 > a/2, situation 2 is 
favored over 1) for any value e > and the depletion 
scenario 1) becomes metastable. However, for a large at- 
traction width, S > a/2, the bridging case is ill-defined 
and is possibly indistinguishable anymore from the sol- 
vated case 2). We note that equivalent scenarios apply 
in mixtures of colloids and nanoparticles i ' 



The findings above may be important to understand 
controversial effects observed in experimental LCST mea- 
surements of polymers and simple peptides. Only if 
polymer-cosolute attractions are weak and not too short- 
ranged, swelling effects should be observable. This may 
be indeed valid for weak hydrophobic or dispersion at- 
tractions. In computer simulations in fact urea swells 
purely hydrophobic polymers! 20 ' 44 In hydrogen bond- 
ing systems, however, the attraction length is short 
~ 0.1 — 0.2 nm and no swelling is possible. Consis- 
tent with that view, urea collapsed PNIPAM polymers 
which has been argued is due bridging by short-ranged 
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FIG. 6: Free energy branches for the situations l)-3) depicted 
in Fig. 5 versus the cosolute-monomer interaction energy e. 
Parameters chosen here are density er 3 p° = 1 and interaction 
width S = a/3. 



H-bonds, strongly supported by experimental meansJ^ 
The same experiments also showed urea-induced swelling 
of hydrophobic peptides. Consistent with our picture it 
was argued that this is due to much weaker, not short- 
ranged attractions.-? 

Looking at the quantitative numbers predicted by the 
present theory they are surprisingly close to the ones 
observed in our simulations. The chain was found to 
be most swollen in the weakly attractive regime e ~ 
0.3 — 0.6 for not too high values of intrapolymer attrac- 
tion. Thus although highly simplified, the effective one- 
component approach seems to capture most of the un- 
derlying physics. 



IV. TWO-COMPONENT FLORY - DE GENNES 
MODEL 



Another perspective to polymer collapse in a highly 
attractive cosolute dispersion could be based on a mean- 
field Flory-de Gennes picture where all interactions are 
described by 2nd and 3rd order virial coefficients. In 
contrast to an effective one-component model as typically 
used in literature, cf. eq. (p}, we now investigate the full 
two-component description. 

Consider a polymer with N monomers and an effective 
bond length b in contact with a reservoir of cosolutes with 
density p®. As usual in Flory theory we take now the end- 
to-end distance R to represent the size of the polymer 
coil. In this view the monomer density profile is just a 
step function with a constant monomer number density 
p m = N/V = 3N/(iirR 3 ) inside the coil with volume V. 
In the spirit of mean-field Flory theory for real chains the 
excess free energy due to interactions can be expressed in 
terms of the virial expansion in the monomer density p m 
and the cosolute density p c within the coil. A Flory-like 



free energy of the system can be then written as 
3R 2 ir 2 Nb 2 



F(R) = 



2Nb 2 



12R 2 



(10) 
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where the first term denotes the elastic free energy of the 
ideal chain, the second term is a correction due to con- 
finement entropy for highly collapsed states^^ and the 
last two terms are the virial corrections up to third or- 
der. The density p c = p c {R) denotes the cosolute density 
inside the polymer which is related to the bulk concen- 
tration p° c by 



p c (R) = p° exp(-^ cxc ) 
with the the excess chemical potential 
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The coupled eqs. (|10)) to (|T2j) provide a free energy ex- 
pression as a function of R, the polymer size for given 
interactions Vy, with i = c,m, as provided by eq. ©, 
polymerization N, bond length 6, and cosolute concen- 
tration p®. To obtain the equilibrium radius R and cor- 
responding monomer and cosolute densities we minimize 
(JTUJ) according to dF/ dR = numerically using an iter- 
ative Newton-Raphson scheme. Due to the relative sim- 
plicity of the equations this takes only seconds on a single, 
ordinary computing processor. 

The R>2 values are explicitly calculated using definition 
(|4|). To connect properly to our computer simulation, 
where e cc = 0.3 and thus B^ c vanishes, we set B^ = 
accordingly. We have also explicitly calculated all values 



B 



ijk 



3 but found that they can be well approximated by a 
constant B l J k — 2a 6 . The latter actually closely corre- 
sponds to the 3rd virial coefficient of hard spheres with 
diameter a. In other words the B% values reflect mostly 
packing effects by hard spheres with size a. 

Results for the polymer size R(e mc ) are plotted in 
Fig. 7 as a function of e mc for e mm = 0.2, 0.3, 0.4, 
and 0.6 (cf. curves in Fig. 2). The radius is scaled by 
Ro = 3.27 nm which is the ideal chain size according 
to the Flory approach. As in the simulation the chain is 
collapsed for small e mc due the addition of repulsive coso- 
lutes. Increasing e mc leads to chain swelling up to a max- 
imum coil size with a corresponding e™^ that increases 
with e mm . This trend reproduces the one found in the 
simulation. Further increase of e mc leads to deswelling 
and reentrant collapse. For the three smaller values of 
£mm the collapse transition is at about e ~ 0.7. This 
value is interestingly in agreement with the value ob- 
tained in the effective one-component approach for the 
crossover from swelling to collapse. This may be acci- 
dental. However, all three approaches show the same 
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FIG. 7: The polymer radius R scaled by its ideal value Ro 
from minimization of the two-component Flory-de Gennes 
free energy eq. (|10|) versus polymer-cosolute interaction 
strength e mc . The plot is for a polymer-polymer interaction 
e = 0.2, 0.3, 0.4, and 0.6 . The polymer shows transitions 
from collapsed states to swollen and collapsed states again , 
in qualitative agreement with the Langevin computer simula- 
tions (cf. Fig. 2). 



qualitative features and numbers in the same ballpark. 
Importantly, also the two-component Flory describes the 
reentrant collapse in strongly attractive cosolute condi- 
tions. Note that the mean-field theory does not provide 
'swelling' as defined by R/Rq > 1. Trends, however, are, 
as discussed above, nicely in accord. 




le I 

mc 

FIG. 8: The mean monomer density p m (solid lines) and coso- 
lute density p c (dashed line) inside the polymer coil versus 
monomer-cosolute attraction e mc and varying e mm = 0.2, 0.3, 
0.4, and 0.6 (see legend). The curves are from the solution of 
the two-component Flory-like model for polymer collapse and 
swelling. For the three smaller e mm values, at the minimum 
of p m the cosolute density p c is comparable to its bulk value 
(horizontal dashed black line). 

The mean monomer and cosolute densities p m (e mc ) 



and p c (e mc ), respectively, are plotted in Fig. 8 for e mm — 
0.2, 0.3, 0.4, and 0.6 as a function of e mc . The monomer 
density p m decreases and increases according to the 
swelling behavior p m oc R~ 3 (cf. Fig. 7). The mean 
cosolute density p c monotonically increases with e mc . For 
small values e mc < 0.4, cosolutes are depleted from the 
polymer coil due to repulsive interactions. At e mc = 0.4, 
the density goes down to one third of the bulk density, 
comparable to the decrease found in the simulations, cf. 
Fig. 4. For the three smaller values of e mm at around 
£mc = 0.4, where the polymer is most swollen, there 
is a plateau in the mean density which almost exactly 
corresponds to bulk density. Remarkably this repro- 
duces the observation in the simulation profiles in Fig. 4, 
where we found a homogenous bulk-like density in the 
most swollen state of the coil. For the largest plotted 
(-mm = 0.6 where polymer compaction is very strong this 
correspondence between the minimum in p m and plateau 
p c is lost for whatever reason. For further increasing 
e mc > 0.5 the cosolute density p c within the polymer fur- 
ther grows for all values of e mm while the chain is collaps- 
ing again. Thus, consistent with the simulations, chain 
collapse driven by highly attractive polymer-cosolute in- 
teractions leads to highly dense cosolute states inside the 
coil. The density of about 14 nm -3 is about twice as high 
as in the simulation, cf. Fig. 4, but the Flory prediction 
of one-order-of-magnitude increase inside the globule is 
qualitatively correct. 
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FIG. 9: Polymer size R{p c ) in dependence of cosolute bulk 
density pi from the two-component Flory-like approach. Data 
are for a fixed monomer-monomer attraction e mm = 0.4 and 
varying monomer-cosolute attractions e mc = 0.1, 0.4, and 1.0. 
Note that only the = 0.4 exhibits nonmonotonic 

behavior. 



Finally, we can use the Flory-like approach to make 
predictions about the dependence of polymer size R on 
the cosolute density reservoir (or bulk) density p®. The 
calculated R(p Q c ) is plotted in Fig. 9 in dependence of 
cosolute bulk density p° c for a fixed monomer- monomer 
attraction e mm = 0.4 and varying monomer-cosolute at- 
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tractions e mc = 0.1, 0.4, and 1.0. The density range cov- 
ers typical denaturant concentration from 1 to 8 nm -3 , 
that is, 1.7 to 13.3 mol/1. In the case e mc = 0.1, R 
decreases monotonically with density. This is expected 
for a depletion-like collapse mechanism, where the effect 
grows with density. In contrast, in the case e mc = 1.0, 
R increases monotonically with density. That indicates 
that the collapse effect due to highly attractive cosolute 
deteriorates with increasing cosolute density. Remark- 
ably, for e mc = 0.4, where swelling is found for weakly 
attractive cosolutes, the size behavior is nonmonotonic. 
This is interesting in the light of LCST measurements 
of poly(NIPAM) and elastin-like peptides in Hofmeister 
saltsi 18 ' 27 Nonmonotonic behavior of the LCST change 
with salt concentration was indeed observed only in the 
swelling scenarios. This happened for the salts Nal and 
NaSCN which seem to feature weakly attractive interac- 
tions with nonpolar peptide groups^ 

V. SUMMARY AND CONCLUDING REMARKS 

The coarse-grained computer simulations performed 
here and elsewhere 2 ^ demonstrate polymer collapse 
to compact globular states at highly attractive coso- 
lute conditions. The polymer response and relative 
swelling/ collapse behavior is most pronounced for poly- 
mers in nearly ideal 0-like solvent conditions. Here the 
polymer can shrink or swell considerably, depending on 
polymer-monomer interaction strength. Collapsed states 
for repulsive and strongly attractive cosolutes are phys- 
ically different as in the latter state the polymer coil is 
rich in cosolutes. For the same polymer size the internal 
cosolute density can differ enormously between 1/2 to 1/3 
of the bulk density (repulsive cosolutes) and 8 times the 
bulk density (highly attractive case), respectively, sig- 
nifying considerably different physical properties of the 
globule-cosolute system. The most swollen polymer state 
is characterized by a bulk like cosolute density inside the 
polymer coil. 

We further demonstrate that the transition to col- 
lapsed states in strongly attractive cosolutes can be 
rationalized by two independent theoretical and semi- 
analytical approaches. Their implications are summa- 
rized in the following. 

The 'effective one-component statistical mechanics 
model' showed that for large polymer cosolute attrac- 
tion bridging interactions (leading to compact 'glued' 
states) are energetically favored over both depletion in- 
teraction (favoring collapse by exclusion) and 'repulsion 
by attraction' situation (leading to swelling). We have 
shown that a critical parameter to distinguish those sce- 
narios is the ratio between interaction width and cosolute 
size 5/o~. Only if polymer-cosolute interactions are not 
too short-ranged swelling effects should be observable. 



This may be valid for hydrophobic or dispersion attrac- 
tions. In simulations indeed urea swells purely hydropho- 
bic polymer o 20 ! 44 or hydrophobic peptides. 19 In hydro- 
gen bonding systems, however, the attraction length is 
short 6/a ~ 1/4 or 1/5 and no swelling is possible. Con- 
sistent with that view, urea compresses PNIPAM poly- 
mers which has been argued is due to H-bonding and 
crosslinking by bivalent binding as demonstrated by ex- 
periments* 1 ^ Methylated urea did not swell PNIPAM, 
neither did urea swell elastin-like peptides. 

The two-component Flory-de Gennes mean-field model 
also predicts polymer collapse at strongly attractive coso- 
lute conditions. The transition from swollen to col- 
lapsed state at a polymer-cosolute attraction energy of 
— In 2 ksT is in good agreement with our computer sim- 
ulations and the effective one-component approach. In 
particular, the Flory model confirmed that cosolute den- 
sities inside the coil can significantly differ between re- 
pulsive cosolutes and strongly attractive ones. It also 
reproduced the simulation finding that the polymer is 
most swollen if internal cosolute density exactly matches 
bulk density for nearly ideal cosolutes. Thus, the struc- 
tural nature of collapsed states can be qualitatively dif- 
ferent while the polymeric state alone, on a first glance, 
may be similar. The Flory approach also predicts a non- 
monotonic dependence of swelling magnitude with coso- 
lute density, which is observed for PNIPAM or elastin in 
Nal and NaSCN i 18 ' 27 A more detailed investigation may 
shed more light on the reason why. 

We note that dense cosolutes in polymer coils may have 
strong impact on internal friction and viscosity of poly- 
meric globules. Consequently, the conformational kinet- 
ics of polymers maybe highly affected by explicit direct- 
binding cosolutes in contrast to isolated polymers. 37 

Our findings may have implications on the interpreta- 
tion of denaturant action on proteins and the nature of 
denatured states. For instance, the structure of molten 
globular proteins denatured by high NaC104 (sodium 
perchlorate) is controversially discussed ^ 29 ' 30 Consid- 
ering our results and the experimentally known strong 
affinity of NaC104 to the protein backbone ) 10 ' 11 ! 28 the 
reason for the difficulties may be a large amount of bridg- 
ing or crosslinking denaturant which probably strongly 
obscures the experimental characterization of structure 
and its classification. 

Finally, similar considerations as above may be im- 
portant in the interpretation of ion-induced collapse of 
polyelectrolytes beyond simple electrostatics £2r— 
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